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Abstract 

In the framework of the relativistic independent quark model the parameters 
of the QCD-motivated static potential and the quark masses are calculated on the 
basis of the 1 meson mass spectra. The value of the confining potential coefficient 
is found to be (0.197 ± 0.005) GeV^ for quark-antiquark interaction independently 
on their flavours. The dependence of the quasi-Coulombic potential strength on 
the interaction distance are consistent with the QCD-motivated behaviour. The qq- 
separations are evaluated and the e~^e~ decay widths are estimated with the help 
of relativistic modification of the Van Royen-Weisskopf formula. 
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The calculation of spectroscopical hadron characteristics within the data precision [1] 
still remains among the unsolved QCD problems. The main difficulty consists in con- 
sideration of the nonperturbative QCD effects, which are not understood yet from first 
principles but only in terms of various QCD-inspired models. Thus, at the present time 
phenomenological hadron models can be used to a certain extent for the data interpreta- 
tion. Among phenomenological hadron models of different kinds ( e.g., Refs. [2-8] ) the 
relativistic model for quasi-independent quarks seems to be one of the most interesting 
and simplest [9]. In our previous papers [10-13] we developed the translation invariant 
version of this model for description of spectroscopical characteristics of both the light and 
heavy vector mesons. We used the Dirac equation with a static QCD-motivated potential 
to describe the quark or antiquark motion in the mean field inside the meson. The model 
does not contradict the first principles of QCD and allows one in a simple fashion to carry 
out numerical calculations of the meson characteristics, such as their masses, the average 
separations between quarks and antiquarks and the e+e" decay widths of vector mesons. 
However, in those papers the effective potential entered in the second order model radial 
equation was chosen for simplicity in a approximate form. In particular this simplification 
lead to the constancy of the quasi-Coulombic potential strength evaluated on the basis of 
1 meson mass spectra for different quark fiavours. 

In the present paper in the framework of the relativistic independent quark model [10- 
13] we use the exact expression for the static interaction effective potential and evaluate 
the potential parameters coupled with the quark masses on the basis of 1 meson mass 
spectra. Then the obtained values are used for calculating the masses of 1 radial 
excitations. The e+e~ decay widths of heavy mesons are estimated with the help of 
relativistic modification of the Van Royen-Weisskopf formula. The gg-separations for 
the light and heavy mesons are evaluated as well. We find an alternation of the strong 
coupling constant ag (which depends on the quark fiavour and is consistent with the QCD 
asymptotic freedom phenomenon) together with the constancy of the confining potential 
coefficient for all quark fiavours. The latter result supports the hypothesis of the fiavor 
independence of confinement potential [14-16]. 

According to the main statement of the independent quark model the hadron is con- 
sidered as a system composed of the non-interacting with each other directly valent con- 
stituents (quarks and antiquarks) moving in some mean field. One supposes that this field 
is a colour singlet confining field, which is produced by the constituents and takes into 
account the effects of creation and annihilation of the sea qq pairs as well. To simplify 
calculations it is assumed that this mean field is spherically symmetric and its motion 
in space is determined by motion of its center. Furthermore, in order to keep out the 
ghost-motion states of the center of mass for the independent constituents, we associate 
the center of mass with the mean field center. We treat the mean field as a quasi-classical 
object possessing some energy, however one cannot evaluate the mean field energy without 
additional assumptions, and this quantity is a phenomenological parameter in the model 
considered. Each of constituents interacting with the spherical mean field gets the state 
with a definite value of its energy. On the phenomenological ground the mass formula for 

jPC 

meson, which is n^^^^Lj state of gg'-system, can be presented in the following form 
[10-13]: 
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Mr,iJ''^)=E{n,.,j) + E{n',,f), (1) 

where nr{n'j.) and are the radial quantum number and the angular moment quantum 
number of the quark (antiquark). E{nr,j) is called the energy spectral function or the 
mass term for the quark and in the framework of our model it has got the following 
phenomenological expression: 

Einr,j) = [A + mT/' + c[l + {^1^^^^'^% (2) 

The first term represents the rclativistic effective energy of the constituent moving in the 
mean field inside the meson, the quark mass m being a model parameter and A being 
found as an eigenvalue of the radial equation of the model. The second term in the 
formula (2) includes a part of energy of the mean field and is purely phenomenological 
one, which gives small corrections to the spectra of different meson families. 

In the equal-time approximation the entire system, the vector qq meson, is described by 
a stationary wave function in the center-of-mass frame, which in the translation-invariant 
model [10-13] can be represented as follows: 

$,(r,r') = 0oVi(r)e^7Mt^cV'(r'), (3) 

where r and r' are the quark and antiquark coordinates in respect to the mean field 
center, 00 is some constant, is a 4-vector of the meson polarization, t] describes a 
meson polarization state and takes the values (1, 0, —1), Uc is the charge conjugation 
operator. Thus, a constituent turns out to be coupled to the other via the common mean 
field which is moving together with their center of mass. 

The wave function il}{r) for any constituent (the quark or antiquark orbital) is a 
solution of the single-particle equation with a static interaction potential. We choose the 
Dirac equation with the QCD-motivated potential V{r ) for description of interaction of the 
quark or antiquark with the external mean field in order to determine the orbitals inside 
the meson. The potential V{r) is spherically symmetrical and consists of the Lorentz 
scalar and vector parts: V{r) — PVo{r) + Vi{r). Hence the equation for a fermionic 
constituent is 

VX + m2^(r) = [(ap) + p{m + Vo) + Vi] ^(r) (4) 

with Vo{r) = ar/2 and Vi(r) = —2as/3r, where the two model parameters a and are 
introduced and have got meanings of the string tension and the strong coupling constant 
at small distances, correspondingly. 

It is well-known that the solutions of Eq. (4) with the total angular momentum j and 
its projection m can be represented as 

where n — r/r. li k — —uj{j + 1/2) where uj is an eigenvalue of the space-parity operator, 
the system of the radial Dirac equations for the fermion in the mean field with the definite 
energy sign and spin projection reads 
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2-Vo-V,-m)f^-^^±^9-9', 



(VA + + V0-V1 + m)g 



r 

[k-l] 



f + /'■ 



(6) 



Using Eqs. (6) one can derive the second order equation for the "large" component 
/(r), and then making a substitution 



r / 1-1/2 

ip{r) — rf{r) Vo(r) — Vi{r) + m + v A + 
one comes on to the model radial equation for (p{r) in the following form: 



(7) 



if + Xcp ^ 



(m + Vof - (VaT^ - V^f + ^^^-^ + 



3(^0 - v;f 



4( VA + m2 -V1 + V0 + m)2 

(v^o" - v;') 



Vi + Vo + m) 2(VA + m2 - 1/1 + 1/0 + m)_ 



(p. 



(8) 



Further on we restrict ourselves evaluating only characteristics of the radial excitations 
of the S-wave 1 mesons because they are described by the simplest version of model 

and supported by the most extensive set of accurate data [1], especially for heavy mesons. 
So, k = 1 and the contribution of the term containing c in Eq. (2) is equal to zero. A 
entering Eq. (2) can be calculated with the help of the S-wave radial model equation: 



if + Xif 



4:asVX + m? f2a 



ar 

T 



[ 3r V 3r 

(m + VA + 7v?^ + 



mar 



3r2 



(9) 



2q;, / n. o\ cr^^ 

— - — h r [m + V A + 1 H — — - 



The right-hand side of Eq. (9) has a singularity at the origin, and when r — > it 
behaves as 3/4r^ — 4Q;g/9r^. Therefore one should keep a, < 3/2 in order to prevent a 
fall at the origin. At the infinity r — > cxo the effective potential behaves as the oscillatory 
one and tends to (7^r^/4. 

In the framework of our model the quantities m, a and as are phenomenological 
parameters and to be determined as a result of calculations of meson mass spectra, as 
well as the e+e~ decay widths, and a comparison of them with experimental data. 

The model equation (9) can be solved only by numerical methods. For calculating its 
eigenvalues we used the computer code algorithm, which is based on the Numerov three 
point recurrent relation for the equation y" — F{r, E)y [17]: 
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Vi-i = 



y,(2 + ^F{n,E)h'^ 



- Vi+i 



1 - 



12 



) 



X 1- 
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• (10) 



For bound states y{r) tends to zero when r > ra- Here rd is a classical radius of the bound 
state, which is determined from the equation U{rci) = E and E is the initial energy of 
the considered level. Then the E value can be determined from the condition: y{0) = 0. 
However, F(r, E) has a singularity at the origin. So, when we calculate y{0) with the 
help of formula (10) we use a regularization procedure for the singular potential, which is 
analogous to the procedure worked out, for instance, in Ref. [18]. 

In order to estimate the precision of the calculation algorithm the well-known test 
with diminished spacing h is used. The results obtained allow us to disregard machine 
calculation errors as compared with systematic errors of the model. The accuracy criterion 
for fitted parameters was the maximum value of the acceptable errors for evaluated hadron 
mass values as compared with typical experimental errors which we choose to be in the 
range of 30 -r- 40 MeV. Thus, the parameter errors written below must be considered as 
our estimation of systematic errors of the model. 

When calculating the values of the model parameters we do not suppose a priori the 
validity of the hypothesis of flavor independence for the confining potential. First of all 
the model parameters m, ag and a for the bb and cc radial excitations of 1 states 
were found. The fit was carried out for each meson family independently. The existing 
experimental data allowed to find these parameters as well as to prove the validity of the 
model. In this manner it was found that the vahic of a is the same for the bb and cc states 
within the systematic errors of the model and equal to (0.197 ±0.005) GeV^. Taking into 
account that there are no well established data for the higher radial excitations of 1 
light mesons, the obtained a value was used for calculating m and for light quarks. 
Note that when doing the calculations we neglect the isotopic mass splitting of the u- and 
d-quarks because it is beyond of the model accuracy. 

All results for the well known experimental data on the 1 mesons, which are com- 
posed of the quark and antiquark with u-, d-, s-, c- or 6-flavours, are consistent within the 
model accuracy with the following values for the coupling ctg ( quasi- Coulombic potential 
strength) and the quark masses: 



Now let us estimate the e+e~ decay widths of the vector mesons with the help of 
solutions of the Dirac equations (6) which were obtained numerically. It is well known 
that in the nonrelativistic potential approach the vector meson e'^e" width is described 
by the Van Royen-Weisskopf formula with first order QCD correction: 



m„,rf = (0.01 ± 0.008) GeV, a^''^ = 0.70 ± 0.15, 

(0.18 ± 0.03) GeV, = 0.50 ± 0.10, 

mc = (1.35 ± 0.05) GeV, = 0.36 ± 0.03, 

mb = (4.70 ± 0.10) GeV, = 0.24 ± 0.02. 




(11) 
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where |\l/(0)p is the square of the meson wave function, which describes the probabihty 
density of meson at the origin. We shall keep a structure of the formula (pi]) in the 
relativistic independent quark model. 

Using Eqs. (3) and (5) the wave function of the meson consisting of independent quark 
and antiquark in the S'-state can be written as 



^(r,r') 



An 



-f{r)f{r') w {(Ter,)w' + g{r)g{r') w (crn)(cre^)(crnV 



(12) 



where we have used a relationship a2 cr= — crcr2 and denoted w' = —a2 w (the asterisk 
denotes complex conjugation). So the probability density of independent quark-antiquark 
system can be easily calculated: 



|<l>^(r,r 



'M2 



167r^ 



{f\r)f{T') + g\T)g\T') + 2/(r)(7(r)/(r')^7(r0(nn') 



- 2f{T)g{r)f\r')g{r')[{e,^){K, n') + (e, n)(e,n')]}. 
And the averaged over initial spin projections probability density is 



(13) 



|$„(r,r' 



levr^ 



[f{r)f\T') + g\T)g\T') + 2/3/(r)^(r)/(rO^(r')(nn')]. (14) 



Now taking into account that r and r' are the coordinates of the quark and the 
antiquark in respect to their center of mass, which coincides with the center of the mean 
field, we can define according to Eq. ([l^ the probability density of meson in the following 
way: 



1 

2^ 



where we denote 



drr f (r) = Ki, drr g (r) = K2- 



Hence the full probability is 



(15) 
(16) 

(17) 



and therefore the normalized propability density which is analogous to the |\E'(r)p is equal 
to (Ki/^(r) + K2g'^{r)) /4tt (k^ + k"^)- This quantity at r = must enter the formula (|TT]) 
instead of the |\E'(0)p. However, in relativistic theory one cannot localized a quark within 
the volume with the radius less than the Compton wave length of the quark Xg = m~^. 
According to this reasoning we use for the estimation of decay width the following formula: 



2^2 



Vn 



1 - 



'1^ 



Taking into account the expression for the normalized probability density, the mean 
squared radius of the meson is 
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p2 

1 



drr 



(19) 



The results obtained for the meson masses and the average radii R^^ are displayed 
in Table 1. We estimate the e^e~ decay widths F^'^ with the help of formula ([T8| ) only 
for the cc- and 66-mesons because there are no satisfactory data for the radial excitations 
of light mesons. The mass values and M^^^ coincide with each other within the 
model accuracy. Taking into account the approximate nature of the formula (^) there is 
a conformity between the calculated decay ratios and the experimental values [1]. How- 
ever, some futher efforts are needed in order to improve the precision for decay widths 
calculations in the framework of model considered. Note that there is a good agreement 
between the i?^^ and the values calculated in the framework of the lattice QCD in the 
valence quark approximation [19]. 

In this paper we considered the relativistic hadron model with independent quarks and 
evaluated the parameters of the qq model potential, as well as the masses, the average 
radii of the 1 mesons, and estimated the e^e~ decay widths. We find that the value 
of cr = (0.197 ± 0.005) GeV^ is the admissible string tension value for u-, d-, s-, c- and 
6-quark flavours within the systematical errors of the model. It is an important result of 
our paper which confirms the flavor independence of confining potential on the 3 ■ 10~^ 
level of accuracy. Thus one may interpret a as the universal parameter of quark- ant iquark 
interaction at large distances and consider the characteristic confinement scales such as the 
confinement length 1^ and the confinement mass rric, which satisfy the following relations: 
mJc = 1, a = ml. 

In addition, the correspondence to the QCD asymptotic freedom phenomenon, or the 
decrease of as when the interaction range diminishes, is confirmed both for the light and 
heavy quarks in the framework of this model. Note that within the large systematical 
errors of the model it is difficult to prove the logarithmic dependence of on Q^, more- 
over, there are considerable nonperturbative contributions for the light quarks (see, e.g. 
Refs. [20, 21]). Nevertheless, there is no contradiction between the as{Q^) behaviour in 
QCD and the diminution of the parameter in this model when coming from the heavy 
to the light quarks. 

An advantage of the presented model is its conceptual clarity and technical simplicity, 
which allows one to calculate in a simple fashion the meson characteristics such as masses, 
gg-separations and e~e^ decay widths with an admissible among phenomenological models 
accuracy. For instance, the values of the parameters mi, and mc evaluated in the framework 
of this model are similar to the values of the pole masses and M^. obtained in Ref . [7] . 
The values of asimf) {i = b, c) are close to the as{Mf) in the MS scheme, whereas the 
average gg-separations for heavy vector mesons are in agreement with those found on the 
basis of the lattice QCD [19]. Moreover, in this paper the calculations of gg-separations 
for light vector mesons are presented, while the same calculations are not possible now in 
the frame of the lattice QCD. In order to increase the model accuracy it is reasonable to 
take into account the relativistic corrections due to contact particle-particle interaction, 
like those which have been treated in Ref. [2]. The account for such kind of corrections 
lies beyond the scope of the independent quark model and will be of our further task in 



Table 1 
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order to improve the model accuracy and to find its limitations. Besides we intend to 
apply this model to the more complicated systems, such as the baryons and the exotic 
multiquark mesons. 
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Table 1. Evaluated masses and e'^e decay widths of the 1 mesons in comparison with the 
data from Ref. [1], and the average gg-separations in comparison with the results 
ofRef. [19]. 



Meson 


[M e V \ 


/I /fth r /I /f ^T/i 
[MeV \ 














P 


1-7/:} r\ n 1 A o 


740 


0.65 






6.77± 0.32 


f 

P 


14d5±25 


1455 


1.12 








n 

P 


1700±20 


1730 


1.32 








1 



1019.413±0.008 


1010 


0.58 






4.43±0.05 


J./ 


1680±50 


1650 


u.yt) 








(P 




zuou 


1.34 










3096.88±0.04 


3060 


0.41 


0.43 


7.32 


5.26 ±0.37 


i,' 


3686.00±0.09 


3650 


0.79 


0.85 


2.03 


2.14 ±0.21 


i," 


4040±10 


4070 


1.06 


1.18 


0.97 


0.75 ±0.15 




4415±6 


4390 


1.30 


1.47 


0.58 


0.47 ±0.10 


T 


9460.37±0.21 


9470 


0.26 


0.24 


1.30 


1.32 ±0.05 




10023.3±0.3 


9990 


0.54 


0.51 


0.37 


0.520 ±0.032 




10355.3±0.5 


10325 


0.75 


0.73 


0.21 






10580±3.5 


10550 


0.94 


0.93 


0.16 


0.248 ±0.031 




10865±8 


10830 


1.10 




0.12 


0.31 ±0.07 




11019±8 


10985 


1.25 




0.08 


0.13 ±0.03 
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